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For n > 4, this follows from the fact that ordinary space time has four dimensions 
and, therefore, more than four vectors are always linearly dependent. The basic 
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tation of the function ASER (x;a) obtained in this way has a kernel consisting of 
a Hankel function multiplied by elementary functions only. This representation 
can, in principle, be used to determine the analyticity domains for all the functions 
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Introduction 


he characterization of a function f(x) with a Fourier transform f(k) = | dre ””f(æ) 

that vanishes, except for positive values of k, is a well-known mathematical problem 
with many applications in theoretical physics. The solution of this problem is a func- 
tion f(x) which is the boundary value of an analytic function of æ, regular in the half 
plane Im (x) > 0. In the theory of quantized fields one is interested in a relativistic 
generalization of this problem, which can be formulated in the following way: What 
are the" properties” of an finvariant= function F(x) depending "on "al four-vectorfæ 
= (77, To, 13, T4 = 1%9) and with a Fourier transform F(k) that vanishes except when the 
four-vector k lies in the forward light cone, Er Minen JE = kå + k2 | ke Kå k? Kå <0 
and k, >0? The case where the function F(x) does not depend on any other four- 
vector except æ is comparatively simple, but considerable complications arise when 
the function F(x) depends also on other four-vectors. 

Questions of this kind are of special interest in connection with the properties 
of vacuum expectation values of products of field operators in different space-time 
points. The mathematical structure of vacuum expectation values of this kind has 
been the subject of some recent investigations?;”%?, In particular it was shown in 
relithatthevacuumlexpectationkvaluefofithetproduetrof m(scalar) helds FAT (æn) TEE 
År (Tx) is the boundary value of an analytic function depending only on the Lorentz 
im vara nt var ab les Er (rn) (er ar ET) ET hs analyticfunetion isfresulartin 
a certain domain W that is obtained if one adds an imaginary vector %, to the coor- 
dinate difference x,—%,…1 and.…lets the m/s vary independently inside the forward 
light cone. This result is a consequence of the following two simple assumptions: 


I. The theory is invariant under Lorentz transformations. 
II. The energy-momentum spectrum of the theory contains vectors only in the 
forward light cone?, 


1 Cf. e. g. H. UMEzaWA and S. KAMEFUCHI, Prog. Theor. Phys. 6, 543 (1951); G. KALLÉN, Helv. Phys. 
Acta 25, 417 (1952); H. LEHMANN, Nuovo Cimento 11, 342 (1954); M. GELL-MANN and F, E. Low, Phys. 
Rev. 95, 1300 (1954). 

2 A. WIGHTMAN, Phys. Rev. 101, 860 (1955); D. Ha1L and Å. WIGHTMAN, Mat. Fys. Medd. Dan. Vid. 
Selsk lenor 5 (1957); 

3 G. KALLÉEN and A. WIGHTMAN, Mat. Fys. Skr. Dan. Vid. Selsk. 1, no. 6 (1958). The last paper is 
referred to as KW below. 

4 For further details, see e. g. the introduction of KW. 
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If one imposes the further condition: 


III. A field operator at a point x commutes with a field operator at a point x” if the 
distance between x and x' is space-like, i. e. if (r—x' >0, 


the domain of analyticity of the analytic functions belonging to the vacuum expectation 
values is, in general, enlarged. 

The special case of three field operators was considered in some detail in KW. 
The regularity domain WM for that case was explicitly computed and shown to be 
bounded by so-called ""analytic hypersurfaces”?, The actual calculation of M was 
made in three independent ways?, one of which consisted in getting an integral re- 
presentation of the most general function fulfilling the postulates I and II above and 
investigating the analytic properties of that representation. The purpose of the present 
paper is to generalize this representation of the product of three field operators to 
a representation of a product of n operators. Our hope is that such a representation 
of the n-fold expectation value will be as useful for the explicit determination of the 
analyticity domain for the corresponding analytic function as turned out to be the 
case with the three-fold vacuum expectation value. However, the actual applications 
along these lines are not dealt with in this paper. We note that the vacuum expectation 
value of n+1 (scalar) field operators can be written as a sum over "intermediate 
states” |z7 in the following way (cf. KW): 


7 z 2 GE 9 7 
ORE ERE EN IDE SS ESS BESKED SEERE et em] OS 


STER 


Here, pW is the energy-momentum vector of the state |2,), 


ÆS N tb ; NES slør 5) 
éÉ ==, Tir SØ SE | and AL£æ—1] | År (er) | 2? FS ERE | År | TOD e (DEL DSE: Ya 


Equation (1) means that, if we introduce the Fourier transform Går "dn (p,...p,) of 


the vacuum expectation value <0 | A; (77)... 4,41(7,41)1 0) in the following way, 


1 Y g i: ØR ) 
<0 Åj (77)... Å 141 (7741) UOvÆ oe] ln | dpi. S 13] Do ES TERE min (pg . Pr) > (2) 


- 2 AV, HALS FR NS", f are Ore D G; N - . 
the function G%: "(P1:- -Pn) IS different from zero only when all vectors p, lie in 


the forward light cone. Because of Lorentz invariance we can therefore write 


Sl . ” i l p 
OS DEDE ONE | FEER (3) 
Kæn 2 |Prx | 


where the function G(px p;) depends only on the scalar products Px Pi = Px Pr Pr, Pr 
and is different from zero only if all po: Hence, we can write 


5 Cf. section IV and appendices I and II of KW. 
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c0 


<O | År (77)... Å141(744+1) 107 = i" | | IT day G(— 247) AX 31 (ér; dx), (eg) 


Ad) GAD) 


2 SOM" ; 
with 4391 (Tx; 7) = SÅ SS Ad, sr ilyg SEE nd (PrP1+ 4x1) IT0(px).  (4a) 
k<l 


Equation (4) explicitly exhibits our vacuum expectation value as a convolution 
integral over a "weight function” G(—ax7) with a "generalized singular function” 
AT); (æx; ax)". The two special cases n = 1 and n = 2 are well known? 


. a ER ERE 
ASE (os ADEL = Fr fæ V, el (Sa) 


am FOT Er AE 
DE Volb= Cunoo 


GX YR 
(VOVE 


(+) Syg 
ÅSE (BES EGE 2250) 


HD (VOV) | My 
| 


Ø (== æn 


OQ=a721t 2075279 1 d29229,  Zy = —X, Ty, (5 c) 
Heksen mee 2 
IX = ål E: ti Bo [ad == dg BEEN . (5d) 


Eq. (5a) demonstrates the fact that the function AS? (æ; a) is the boundary value of 
an analytic function regular for all values of z not on the positive real axes. In a similar 
way, Eq. (5b) shows that 4” (7%; ax7) is the boundary value of an analytic function 
of z,; regular and exponentially decreasing for large values of |2Z,,| for all points 
not fulfilling the relation 


Q=YR=0; 0<9=<00. (6) 


It was an exploration of Eq. (6) which led to an explicit determination of the 
boundary of the domain Pt in KW. We are particularly interested in obtaining the 
generalization of (6) for the case n>2. It should perhaps be mentioned that a re- 
presentation of the kind we have in mind generates the most general function ful- 
filling assumptions I and II above, but that the restrictions imposed by assumption 
III are not taken care of. However, it is a separate problem to investigate the en- 
largement of the analyticity domain that follows from assumption III. (For an illustra- 
tion of this statement we refer to sections V, VII, and VIII of KW.) 

The two functions (5) can be calculated by straightforward methods, but a 
frontal attack on the higher functions leads to considerable formal complication. 
Fortunately, there exist relations connecting functions with different values of n as 


6 Functions of this kind were first introduced by A. WIGHTMAN and D. HALL, Phys. Rev. 99, 674 (1955). 
rens rap pendisklsol  ESWE 
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well as invariance properties of each function 4, for a given value of n. The ex- 
ploration of these things will allow us to simplify our calculations to a certain extent. 
The first three sections are devoted to a discussion of these features. 


I. Reduction of the Function 4), to 4/0) when n>5 


Ås ordinary space time does not have more than four dimensions, the vectors 
Px appearing in the definition (4a) cannot be linearly independent when n > 5. This 
allows us to perform some of the integrations in (4a) in a comparatively simple way. 
If we suppose that the matrix || ax, || has rank four, we can choose a set of four 
linearly independent vectors among the px''. It is then only a question of labelling 
to call these vectors p;...p4. With the aid of the d-functions we can then express the 


other vectors ps...Px, as linear combinations of p;:..p4 in the following way: 


n 
SE 
Px 2 Y6A PA fork == DE (7) 


The coefficients &,4 are determined from the conditions 


[old = Og oe Je Do HER. (8) 
This yields 
4 
an ÅR 
BE DEG (9) 
AE 


with 


y= | Fr Fra Ch Con E) (9a) 
dz4 dog dg4 Qd44 


and with 427, being the cofactor of az, in D. Note that the indices 4 and 4" only take 
values from 1 to 4 in this definition. 


After these preliminaries we can now evaluate integrals of the type 
” 


4 
T= der 116 (PrP1 + ax) F(Ppx), (k> 5), (10) 


where F(px) is an arbitrary function of the vector Pr: We get 


1 4 


1 
I= Rk S ax2422 P2' fi) kl 
VED TAP FE (2 
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By repeated application of Eqs. (10) and (11) we find 


2 bd BOLTE R 
Jos NES aDE .dppe'"" HG LSK SUK) Px, + dx k,) 
e . helse k,<k,= 
i Å 1 4 4 (1%) 
II aL -—— )' 411.94247 7a J (5 
re VED ESE ek kk, SE RS LEA GA RAE AAR ENS AAK 


ir 4 
RØRET DET TEA YET AGED 
SD yen KAZAN P2 EB, 


The arguments of the d-functions can be somewhat simplified with the aid of 
the relation 


4 
re: dy År 3,= D 0,3, (83) 
Fl 
if both x and Å are < 4. In this way we get 


=> sy Åan» 
di 5 Dy' 
IR= (FFD) ODESSA =1 Så. rr (cu D- D — adm (Ng (14) 


"kK<Sk'=5 


We now return to the integral (4a) and first remark that we can replace the product 
n 

IIO (px) by 06(p7)II0 (a,,). This follows from the simple observation that the scalar 
k=2 


product of two timelike vectors is negative if they both lie in the same light cone 
and positive if they lie in opposite cones. Using this and (14), we can write the integral 
(4a) in the following way: 


(7 VR 12, P% 4 
AL BEC; 2) = 110 (44) (2 an] del aser TT 0 (Pr Pa) Fa) 4 (p7) 
= ES 
2 iD, 7, | 
x UD le are + OD TT Pr, Pr, + dx x,) | (15) 
. . Make, = 
ÆDES OY n 2 (+) | 
2 + 
E ES IEEÅ D) z Sy mme E 0 Adan dN' ke, VA (vs ) 11044), 
where 
4 AR 
DE: 2: SS IGEN) 


Note that the indices x and x" in JYSK (ur tese konlygtakervaluesktromElktor RE ge) 


shows explicitly how the general function MSP GeKaEcansbelexpresse dase Nara) 
9% 
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multiplied by certain factors. The d-functions appearing in (15) are an expression of 

the fact that the vectors p;...p, are linearly dependent on the vectors p;...p4. 
This result tells us that apart from the two functions exhibited in Eqs. (5) we 

have to discuss only two more functions, viz. 40? and A9?, which somewhat simplifies 


our task. 


II. Connection Between the Functions 4(, and 4/7 when n<4 


In the preceding section we have demonstrated how one can express every 
function 407 with n>5 in terms of the function 4. The mathematical reason for 
this reduction was the fact that space time has four dimensions and therefore more 
than four vectors must be linearly dependent. We now want to remark that one can 
also express 40 with n<5 in terms of 46 with one or more of the vectors æ put 

equal to zero. In fact, if we put e.g. xx, = 0, the result of the p4 integration can only 

depend on the masses ax. This Es from simple considerations of Lorentz in- 
variance. Therefore, the function 47 with xx, = 0 must be equal to the function 
AZ of x,, Yy and æg multiplied by a certain function of the masses. To see this in 
detail we specialize the definition (4a) by putting xx, = 0 in it and obtain 


5] 


(-iy' ? LSD. Yy É 
LE. red de des TT (pe Pr 4 OPD 0 (48) 1, (6) 
FO CAT) k<1<4 k=2 
where 
. 3 1 
Mn = | dr 9 (På + 440) 0 (P4Px + dak) EDER (16a) 


with the determinant D defined by Eq. (9a). We assume explicitly that this deter- 
minant does not vanish. Otherwise, the function AP does not exist. 
We thus obtain the following connection between the functions 4/7 and AP: 


Malle SES TEÅ D) 0 (a47) 44” (æ; a). (17) 


The function dne in (17) depends only on the vectors æ;, xx and æz and on the masses 
ax With k and I both <4. 
Relation (17) enables us to find the function A/? directly from AD, Therefore, 
the only function to be discussed before we know AG, for all values of n is ARD. 
However, it will be useful to consider also the special case of A/ when two 


of the” vectors x;...æ4 vanish, erg. x; ='x= 0 From (4 a)kwetthen find 
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BRED: Bam, > 
SANS kt EN SRENV | | dp; dps e 0 (Pyt 411) 0 (På + 422) 0 (P1 På + 442) 
ROCA) (18) 
; ' 
x 0Q(p3)IT0(a1 4) 134, 
k=2 
where 
== | | dps dp4 0 (pz + 33) å (på + da): ms (Pa Pr t dak) ig. (P3 Pr + 431) 
iÉ 
DEN (=D) | dps å (p3 + 433) 6 (P1 P3 + 413) & (P2 P3 + 423) -(18a) 


3 7 DN ODA: 


Va2, — d1]1 d22 == 


with 
rt, nor ig) 
De = | d12 22 (23 |. (18b) 


d13 d23 (33 


The integral 1,, does not exist if a%, < 444 422. When the determinant D does not vanish, 
it is always possible to relabel the vectors p, in such a way that the determinant D, 
and the expression a7,— 414 422 are different from zero. For simplicity, we assume 
this relabelling to be made in (18). 

From Eqs. (18) we now obtain the following connection between the functions 
Me kande LE 


HED) DYG (a14) 0 (13) AD (11,75; 411, d29, 412) - 


E: =74 =0 (2) Var ar a; dg NED 


In fact, Eq. (19) will turn out to be a very useful tool in our later discussion 
of the function AP. 

For completeness, we want to mention that one can also get a relation between 
the functions 4/7 and 44” similar to Eq. (19). However, this relation will not be 
used in the following, and we do not want to give it explicitly. 

We have now established all connections needed to determine the general 
function i pydk for an arbitrary value of n from the function JA KERDE The following dis- 
cussion is therefore devoted entirely to this function. 


(19) 
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III. Transformation of the Function AG ) to a Standard Form 


We now turn to a discussion of the invariance property of the function ASE 
mentioned at the end of the introduction. For this purpose we remark that we can 
make a transformation of the variables px in the definition (4a) with a non-singular 
real matrix Å,; in the following way: 


4 
Ur= XM 7093 De AT EAOR (20) 


væl 


Note that this is not a Lorentz transformation where the components of each vector 
are transformed among themselves, but a transformation among the vectors where 
all components are transformed in the same way. The scalar products of the vectors 
qx åre given by 


dx dt dy = 0, (21a) 
with 
i 4 
dk ED dy År. (21b) 


If we introduce a matrix notation and write A for the matrix a4;, Eq. (21b) 
can be written in the condensed form 


ÅSE ASAP AGE: ere) 


By elementary considerations one finds the following two formulae: 


Udg sold = (GE | ANE køre. 5 sols (22 a) 
4 i 4 
IT 6 (Pr P, + ag) = (Det | A |)? IT 8 (47 4,4 a47)- (22b) 
pels k<1=1 


Therefore, we can write the integral ÅD as 


4 
12 Y,dr 4 , 


IØ | RER 4 
5 & TE ) sås | dg, .… dg; e IT 6 (dx a+ 4%) 11 0 (ax) 0 (CF (DEDE (23) 
k<l1l=1 k=2 


with 


Ur D Xy GE (23a) 
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Apart from the factor Det | Å | in front and the behaviour of the step functions 
0 (Px), the function AC? is therefore invariant under the combined transformations 
(20) and (23a). 

We now specialize the general matrix Ax, to the following: 


Kkæl 
dx = Pk, rt Pi: (24) 


1=1 


Provided the vectors px are labelled as indicated in connection with equation 
(18), we can determine the constants c,, in (24) in such a way that all the non-diagonal 
products qx 9; = dx; Vanish. A simple consideration shows that the diagonal masses 
ark for that case are given by 


mn 475 (25a) 
É aa a 
me 11 422 » (255) 
ad11 
D 
eee (25c) 


SEE 25d 
BE De (25 d) 
ASEDER| PAN ETS forftheltransformationf (Q4YEWwertind fromt (23) 
(+ 1 Aunas å ag 
Ås == OD DE (LÆKCNSAVE ere … dq4e BEN (97 + DD) OL (1): (26) 


The vector gy is the same as the vector p; and therefore timelike. As all the other 
vectors qx are orthogonal to the vector qg;, it follows that the masses Rim (ØSNEFOD 
k+ 1 must be negative. We then make a new transformation changing all the vectors 
gx by constant positive factors so as to make the absolute value of their squares equal 
iomones bhiskyields 


8 In a space with a Euclidean metric instead of the Lorentz metric, integrals corresponding to our 
AS but defined without the step functions O(P,) are intimately related to the generalized Bessel functions 
studied by S. BocaNER: Med. Lunds Univ. Mat. Sem. Suppl. (1952), p. 12. A symmetry property of these 
Bessel functions corresponding to the invariance of our functions under the transformation (20) and (23 a) 
has also been mentioned by BOCHNER in his paper. For the functions ASB (x;a) studied here, the same re- 
sult was stated and proved in the thesis (unpublished) of D. Ha1r (Princeton 1956). This paper also contains 
the remark that all these functions with n >4 can be reduced to what is here called AS (æx;a). We are in- 
debted to Professor A. WIGHTMAN for making this thesis available to us. 
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1 . 4 
Aj = =D 0(-D) 0 (Dg) å (a72— 411 de KELD (15) 


E (27) 
FE ! SN GEDE GENT (7 ADR 
same Å da da" "ala + 10 (49 718 (d—1) 71 8(42 4). 


where the vectors y, are given by 
ERE ER 1 
Ur > VI! xx | IT) (59 (27a) 
En 


with the matrix Å of Eq. (24). 

The transformation (24) is not the only one that will make all the non-diagonal 
masses equal to zero. In fact, we can make a further transformation among the q4's in 
(27) with a new real matrix A” that leaves the quadratic form q% — q2 — q3 — gå invariant 
and with determinant + 1. With such a transformation the matrix yx, y, Can in many 
cases be made diagonal. We are then left with the integral 


1 1 æ 
AD = Er VED) D DE (7 ae.) I Ca) 28 
5 CD Y-D (Dg) & (a39— 411 22) få (a+) (28a) 
Q BYTE? n= 4 
EN Ga FDR NEO BES DDD Cpepe (28b) 
& = SL 


In the expression (28b) the vectors %x, are all orthogonal to each other. The squares 
of these vectors are determined by the eigenvalues 6;...04 of the matrix Y defined by 


VRE EU SEE KE ERE STE | 
gs URE. JESS (29) 
VE ir. | 
We have 
æ7+0,=0; RED, KEE (30) 


According to (27a), the matrix Y can be written 
FN SE (31) 


The elements Xz, of the matrix X in (31) are given by the scalar products Cy Cy 


1 


of the original vectors x, in (4a). |/—A" is a diagonal matrix with matrix elements 


given by i|/a14, for the first element and by væ Axx for k=1. The eigenvalues of Y 
in (31) are the same as the eigenvalues of the matrix 
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AE 4 SUM READ TT (22) 


where Å is the matrix of the original masses ar; according to Eq. (21c). From (32) 
it follows immediately that the quantities 6, in (30) are the eigenvalues of the matrix 
product — XA. This result allows us to compute the o's directly from the quantities 
appearing in Eq. (4a) without going through all the transformations in detail. 

The condition that the matrix A' transforming (27) into (28) is a real matrix 
is the condition that all the eigenvalues 6, are real. For given matrices X and ÅA this 
may or may not be the case. Instead of trying to discuss explicitly the case where AA” 
is not real, we use the result of ref. 2 that the function 4? is the boundary value of 
an analytic function of the scalar products x, 7æ,. For a given matrix ÅA, we choose such 
matrices X that the quantities 6, are all real and compute the function 42” for that 
case. The value of AG? for some other matrix X can then be obtained from our result 
by analytic continuation. 


IV. Calculation of the Integral (28b) 


The discussion of the previous sections has reduced the computation of all the 
functions AM), to the evaluation of the integral (28b) with vectors x;...æx4 that are 
orthogonal to each other. Further, we know from the result of ref. 2 that this integral 
is the boundary value of an analytic function of the squares æxz. Therefore, it is suf- 
ficient to compute it for, say, æ; timelike with positive time component and hence 
the other æ, spacelike. To make the integral convergent we further assume that %;, 
has a negative imaginary part, which is not necessarily infinitesimal. It is then a 
question of labelling to suppose |æ,|>|x3|]> |x |> 0 and to choose the coordinate 
axes in such a way that %x, lies along the positive æ-direction, xz along the positive 
y-direction, and x4 along the positive z-direction. In this way we get 


va | bål | dn e P10710t IP, Ty yt IPsy sy t ID 4, Z4 5 (D. 1) 0 Cp7) | 
d å 7 . 2 » (33) 
x IT 6(pr,—1) IT 6 (PP). | 
k<l1l=2 


k=2 

The actual order in which the following integrations are performed is a question 
of convenience. We have found it convenient to start with the pz and p4 integrations 
and the integrations over those components over p; and ps» that are orthogonal to 
the plane spanned by æ; and x,. This can also be so formulated that we first compute 
the integral 


UC A 4 
Iy= W dp; dp> dps dp4 6 (p3 DUE GORE Dog oe (34) 
Fz 7 b=2 k<l1l=2 
(y-2) (bx) 
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Note that the integral I, does not contain any exponential functions and, there- 
fore, can be computed by elementary means. Ås is shown in Appendix I, the result 


of this computation is 


9 
ti 


2 


I = 


så 
2 


d ((q1 42)? + (43 44)? — 42 42 — 43 4å) å (43 + q4 — 92 + 91) | (35) 
x 0(A) 0 (92 (g2 + 47) — 2 (91 42)» 

2=(91927—9192>0, (35a) 

JN CD DK CEED CIK DE (35 b) 


In Eq. (35), the vectors gx, are two-dimensional vectors built up of those com- 
ponents of the vectors px that have not been integrated over in (34). To be more 
precise, gq; and q2 have components in the æ and t directions, while g3 and qg4 have 
components along the y and z directions. 

We next perform the integrations over q3 and qg4 and get, according to Appendix II, 


(u» 2) 
ked ner eee SE re Een 
== ln RE i 2V3) SO E 23) 
AE 2 5 (36) 


(fer te gg q 2V3) (2 g GÅ | Vi)! | 
| 


rd) 


x Ø(A) 0 (47 (43+ 43) —2 (41 4277) 9(42—43—2 1/4) 0 (42— 47 + 2V 2). 


According to Appendix II, the integrations over gq, and qg2 can be arranged in 
such a way that the integral I in (33) reads 


2 9 


sd 


f een, NEUE VRE DØRE ap 
g FEY ØER By LÆR By Az By " 42 
ale] Ea En SS RR RR 


The function F(1,, 7,) in (36a) is a certain function of r, and r,, but independent 
of the x,'s. This function is given explicitly in Appendix II, but the details of it are 
not important for the following discussion. 


IL er= BSN, eg 
fe Å Jr Klee (rd) HEN |- SE "ég ri) ES (tot: er | 
FÆDIDNN 
(36a) 
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We now make use of the formula? 


1 IT 
HDD D dvd G), (37) 


5=1/2 42 £2272 008 p, (37a) 


where the sign of the root in (37a) is defined so as to make 3— z for |2|— &. Eq. 
(37) is valid as long as |z2|>|Z|, in which case 3 never passes through the origin 
when y goes from 0 to ax. Therefore, it is important in our application to have the ab- 
solute value of the argument of the Hankel function larger than the absolute value 
of the argument of the Bessel function. One way to achieve this is to make the imaginary 
part of 7, sufficiently large and to compute the integral for that case. We can then 
again use the analyticity properties established in ref. 2 to continue our result to 
arbitrary values of the vectors involved. 
With the aid of (37) we can write (36a) as 


1-50 +19), (38a) 
re 07 7 0 i ; 
ANE gl É a ad TT ED il. Ø 
Kik ie dy; dw> ole TE (Te RENDE ek i 372 32 Å .(38b) 
9 0 


2 [; É 2 and ” ge 
sie = ne GE) (TEE TE) (CE COSTER DEO) 


The term IØ) in (38a) is obtained from IY if we replace e.g. x4, by —X4,. 
By putting x,,= x4,= 0 in (36a), and using the result of Eq. (19) together with 
Bast )Fand 8) ver setter formula 


sg X 
FIS == C19 FE 265, Ée DE To 
w 0 År 9 
t (39) 


iz? 


ig É ik i 5 
(Hf . (Te DR) HG (57 or) 


10 Toy 


Eq. (39) is correct when the imaginary parts of the arguments of the Hankel 
functions are positive. Otherwise, the integral on the left-hand side is not convergent. 
Eq. (39) allows us to perform the two r integrations in (38b) explicitly, yielding the 


result 
ant 
so HØ (- ) JEIND HE 
Id TEN? dw> 2 E. = 52) (40) 
we, BL 1820) 


9 Cf. e. g. G. N. WATSON: Theory of Bessel Functions, second edition, (Cambridge 1944). Eq. (37) 
is a special case of formula (16) on page 367 of this reference. 
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Å simplé algebraic calculation shows that, if we introduce the argument of each Hankel 
function as variable of integration instead of one of the angles w, the integration over 
the other angle p can be made by elementary means and yields 


t2 ta 
KARET EKG ÆRE DAGE 
10 8 ird]] HE) | : 0 (1) i (41) 
METRO STORT AKA ET 
PODET ENE FE) EP DEERE NEDE (41a) 
lt T10 7 Tag" TyyT Tags 15=> —T10— Tag tT3y + TA 2 1 
12 19 l Toy Tsyt Ty; 4 — — 0 FT Tag FT3Yy TT Tyg- | 


For the particular case of the vectors xx, which we have investigated, all the 
numbers fx, have the same positive imaginary part, while the real parts fulfil the 
inequalities 

RAVE kerse=RErE- REE (41c) 


Finally, all the square roots appearing in (41) are defined to have positive 
imaginary parts along the path of integration. When w goes through real values, the 
two paths of integration in (41) are two arcs of hyperbolas, as shown in fig. 1. 

We can now redefine the square roots by introducing cuts in the complex t-plane 
between the: pairs of points (8,8); (17 ==); (571) Tand (FLE) Sandidefinme 


VP, (t) to approach ?? for large |t]|. 


Fig. 1. Paths of integration in Eq. (41). 
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As shown in Appendix III, this new definition of the square roots implies that VP, (t) 
is the same as in Eq. (41), while the sign in front of VP, (t) has to be changed. Further, 
one has to replace V-P, (DEby LÆRE (H) for k=1 as well as for k=2. This gives 


SKOV LOE tdt HØ (t) | 
DEG ON ANE 


ta 


(42) 
(2 NE URET ED) 
red kr (03) 
EDER D DB (UD: (42a) 


/E(t) in Eq. (42) is defined with the aid of the cuts mentioned earlier and the 
condition that it approaches t" for large values of |t|. A straightforward algebraic 
calculation allows us to express the polynomial £(t) in terms of the squares of the 
MEecLOorTSs Xr, 1. e€., in terms of the eigenvalues 0, mentioned in Eq. (30). The result of 
this calculation is 


NE (FO EREVT FPS 8 VT EP (250): (43a) 
R=4 (071094 01034 0704 409034 09044 0304) = 2 (HF —1), (43c) 
BL mr Å 
S= 64 (07 02074 07020440; 0304 4 090304) =—7(11—31, 1,421), (43d) 
8 9 F; & = 
RES [AO EET (43 f) 


where the matrices AÅ and X are defined in connection with Eq. (32). Note that T is 
the product of the two determinants of the matrices &Æ and Å (apart from a numer- 
ical factor), $ is a sum of products of 3x 3 subdeterminants of XÆ and 4, while R 
is a sum of products of 2 x 2 subdeterminants from the same matrices. This means 
that T=0 if the vectors xx lie in a three-plane, S= T=0 if they lie in a two-plane 
and S=T=R=0 if they are all collinear. 

We now finish with the remark that I in (38a) is obtained from IW in (42) 
if we replace Br in (43a) by yo It follows that 
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ODS fra ue0] lø D) 0 (D,) 


pr 
OP VED JE VE E, 
4 


x0 (a32 miOflat Tea) 9 (212): 


(44) 


eee 70Y REV SP SES Y TE 0) (44a) 
where the determinants D and D, are defined in Eqs. (9a) and (18b). 

The result (44) has been derived under somewhat special assumptions about 
the matrix X, but it follows from the analyticity properties mentioned several times 
earlier that it is valid for all X. Therefore, it yields the desired result. In particular, 
the two equations 


are the desired generalization of Eq. (6) for the function AG), with 4>n>2. Note 
that (45) is reduced to (6) when S=T=0. 


V. A Simplified Version of the Result of Section I 


The result of Section I allows us to express all functions AG), with n>4 in terms 


of AG? with the aid of the formula (15). An explicit application of this formula is 
rather involved as one has to compute the vectors y, with the aid of (15a). In this 
section, we want to give an explicit expression for the quantities Q, R, $, and T in 
terms of the original vectors æx when n>4. To this purpose we write Eq. (15a) as 


Uke EVER (46) 
wherelgkandkærarekte eematrices Stisen (ære or) Fame 
is a 4x(n—4) matrix defined by 
1 4 
Mer 2 Ana dar: DE DRÆN ES RE ng (47) 


For convenience, we split the matrix ÅA in the following way: 


Å 9 
1 


År 5 
ABN (28) 
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where A, is a 4x4 matrix with matrix element 4,4 (2,4=1...4), A; is a 4 x (1—4) 
matrix kwithfmatrixtelements "a77 (+ == EAN RES En etc In this matrix notation 
kekcan kvadet Rg ES as 

AM=A4. (49) 


From the å-functions in (15) it also follows that 
AR M == Ås. (50) 
Eqs. (49) and (50) can be combined to yield 
MP AM=A2. (51) 


We now consider the tensor Fiv with two vector indices u and v in ordinary space 
time defined by 


4 
E PREDY 41 (Y2)v = by Å 0 Yv: (52) 


L47) 


With the aid of (46), (49), and (51) we get 


SEERE ØM EE HYESD) 
mn lene 97 ze Ål HE SR RH RS Eg ; (53) 
n 
S RB (Xx)v dx (77)4: 
ren 
The important quantity in the formulae (43) is the matrix GW defined by 
4 4 
GQ=5X S (Yddy Ure Iye 4 EØS (54) 


and the quantities Q...7' are obtained from traces of various powers of GY, However, 
fromt 62 )Fandt (63) Holowskergs 


4 
Sp KG AE DE mu TT = Sp [G97 (55) 
uM=1 
with 
4 n 
GØYE Eg 2 (rd (Cr ut dr = OGRE Tese == I orde (55 a) 


In a similar way, one finds 


== SDP YT So ES sd (56) 
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Eq. (56) allows us to compute the quantities /, and, hence, Q...T directly 
from the given matrices X and Å without going through the intermediate steps of com- 
puting the Y. This is the desired result. Furthermore, we remark that Eq. (56) is 
formally identical with (43f). Therefore we can consider Eqs. (43) as the general 
definition of Q...T for all values of n. 


Appendix I 


Consider the integral I, defined in Eq. (34) 


mm UL 4 
IT = NW du, du, UW du, du, 0 (u3 +1 + 0 (uZ= 1 + q7).0 (U, u24 q192) 
"X ze 


T) 


(w-z) ( (AR 


Xx 0 (Ug U4 + 9394) då (Uj 93 + Ug q1) å (U1 q4 + U4 41) 0 (U5 93 + U3 92) då (Us q4 + U4 q2) - 


In Eq. (A.1), we have written px, = u%x+ 95, Where gx, are those components of 
the vectors p, over which we do not integrate (cf. the remark after Eq. (35)). 

With the aid of a transformation similar to that in (24) we first introduce 
unit vectors eé, orthogonal to each other in the following way: 


q1 93 

SE NES 58 es = Sa. 

| —q1 y 93 

2 2 q 2 / 2 CAR 2) 

q2V/ — 43 7 — Å RR qAY VEL 3 
3 ad V 93 
2 SEE CA gg , 

V742 V 734 


are == BØGER ED 289 2 . ; 
where 412 = (91 92)" — 97 42 and 224 = 43 q4— (93 94)”. If we simultaneously introduce new 
variables of integration according to 


Wu Ulg 
DD] 5z Så Vg = =3 
V- 4 Vq2 
2 ÅZ3 
Us) — 977 i DB i U4 | Q2— ale $ ( ) 
SR V-qi V 43 
2 == ; V4 — 
V2;2 V 434 


the integral I, can be written as 
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— r— U OU 4 É 
1, =V2V SAN dvy dvy Å dvs dv, IT 8 (Vi + Qx) d (Vy Vo + 112) d (V3 Va + 134) | 
(2) (æt) meg | (A.4) 


xO(V;eg3+UV3e;)0(V5,634V362) 0 (V, 8440461) O(VvV5 044 V4 €7), 


with 
1+ q7 (9192) =(97)" 931 (439437 + (43) 
SE an ER 27 —; Er; OUr= SEG] 2 (A.4a) 
q1 — q1 412 93 93 "34 
and 
RE DU CALAB) 


12 == »… 1]34 Fr 
VA dg V 434 73 


As the space of the vectors e, and e, has a positive definite metric, we can perform 
an ordinary rotation in this space and make the vectors vg and v4 orthogonal to each 
other. (A corresponding real rotation in the space of e, and es, is possible sometimes, 
but not always.) In this way, we get a formula similar to Eq. (A. 4), except that 134 = 0 

re (4 [4 . 
and Q3 and 94 are replaced by two other quantities 93 and 04 determined from 


2 2 
' ' « 9379 É 
03+04=03+704= 2 A z GO) 
34 
vø. g  (1—43) (1 — 92) — (93 44) 
Q3 Q4 = Q3 Q4 1184 Er 1 . (AE) 
34 


If we now introduce a coordinate system with the vector e; along the time axis, 
the vector e, along the æ-axis, etc., the last four delta functions in Eq. (A.4) read 


0 (Vi y —V39)  (V2, + V34) Å (UV, — V49) O (V5, + VA 7): FALD) 


These four delta functions permit us to make, e. g., the v4 and vs» integrations trivially, 
yielding 
Io = 47 3 Å dvs dv, 0 (vs F03)0 (vi + 0%) 0 (V3V4) | 
(æt) —— (UA) 


BR > RENE 2 
x 0 (V30+ V40 +01) 0 (V37 + V4x + Q2) O (— V39 V34 — V40 Vag FN 12) - 
The last delta function in (A.6) can now be written as 


0 (1127 V30V37— V40 VA) = 21742] O(Ni2— (V30V37t V40 VA) ) I (N 72 (V30 V3 gt V40 VA 7)) | 
(GA 7D) 
6 2 ! £ £ 
DN ro (Go (EEK OD EEN DE). | 
4x+ 
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In a similar way, we write 


å(v2,+v2,4+02) =0 (027017 03— 24) (A. 8) 
This gives 


9 | VAD V4z4 å (fe = (0; 403) (07 £.04)) d (0203 03 04) | 
i (A. 9) 
x Å dvs dv, å (v3 + Q3) å (v4 + 94) å (Vs V4) då (V39+ vo t Q1) 9 (132 (U30 377 V40 VA 7)): | 


The remaining integrations in (A.9) can now be made by straightforward methods, 
yielding 


Øg 


Va, 


y 3 Q4) (2 Gin (CE 03) (0; + Q4)) MC= 02 Qa) UNDER ØNE, (10) 
— Q3 04 


(0270170 


Q . . / nn , 
If we introduce the expressions (A.4) for Q;,...04 and 142, we get after some 
simple manipulations 


( DU N ere 2 2 2 É « 
I, = DE O(A2—434)0(43+44— 42 + 43) 9 (4) 0 (47 (47 + 42) —2 (41 42)),  (A.11) 


Å= (92=1) (GE DS (49) GE) 


Eq. (A.11) is identical with Eq. (35) in the main text. 


Appendix II 


With the aid of the result of Appendix I, the integral I, in Eq. (36) can be evaluated 
by straightforward techniques. We write the definition of this integral in the following 


way: 
n 


Ch 9 9 9 ? g 
TI = VA 0(A) 0 (47 (42+ 497) 2 (91 42P)I, CARIP) 
5 N dgg dg; ePevtartidis%4r 5 (g3 + gg —07) 0 (93 4 (9394) 4): (A.12a) 
(GREN 
== GE AE (EO O (A.12b) 


The argument of the last delta function is simplified if we introduce a new vector q4 
with the same length as q4, but orthogonal to it, i. e. 
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, 
Jay 7 daz> 


; (A. 13) 
Jazz = — day: 
With this definition we get 
12 
Gerda (48138) 
(44 43)? = 43 92— (9394) = 4, (A.13b) 


and write the integral Jf as 


1 er fg 55 RE Enø pe 5 / An £ / HF 8 
= E /AN dq3 dq4 e"13YT3y T I04, Ta 2 9) (9 +- GE == Ø) (lå (93 Ga 4) ES) (GE VE VÆN in (A. 14) 


If we now introduce the two vectors pz = 43494 and p4 = 43—94 as. variables 
of integration, we get 


1 OFC ID KON LENE Een 
Fr 3y 9 2D 4 , 2 
TR) des dps? gj 3 (re 022 VANG (pE ø eN 


2 


| 
«do(2er 2 zl TAB (rr 7 SE r2yi)/ Ju tt EDP) | 


SO (02 VA) 0072 VA), 


S FØR > z == m Cay tt X47 5 
Ed RHØR ALA sy. STÅ RER 
2 (A. 


where 9 (æ) is' the ordinary Bessel function of order zero. Collecting Eqs. (412) 
anna SS) vere RETTE GO) Fink hermamnetext 

Essentially the same technique can be used for the integrations over gq; and q% 
in (36a). ÅA typical term to be computed is 


1, =Vda, dg, em? ASSR me 000 - ad Fo VG (G= å—2|/4) 


» (A.16) 


æ. sr, Fo GE tg, RL ES Fo SES REE) PA 2 6 ER 
ZA! Te NER arv) mV g2æd? STATEN SETE 2): 


In this case, we choose py = 41492 and po = q1—9; as variables of integration, 
where q, is defined by 
[4 
d2x — d20> | 


få 
q2o = T d27" 
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This yields 5 
(k == =0B- (A.18a) 


(427) (919) 91 72 74; (A. 18b) 


PN ØER To 72, 2 ny Q Å) 
ir P1- P2 e Q 2) Q9(—D2)0 EG 
| 4 ) VA (-p1) 9(=p2) 0 (P1 + P2) 


G) | (A. 19) 


SER og JC SEES DØD: 
Lr 2LZ É 5' 
sJyf tu ter V=R)vo( =V—)e( [Dyd ve dl 


The two vectors p; and ps are both timelike with their sum in the future light 
cone. The scalar product of two timelike vectors p; and ps, is either smaller than 


—|/p? på when p, and ps; lie in the same light cone or bigger than |/p? på when they 
lie in opposite light cones. The last 0-function in (A.19) tells us that 


På på V 

il 2 2 

jog” 2 =VAPS PS (A220) 
sad Bi P> AE 


Therefore, p; and pz both lie in the same light cone. As the sum of p; and pz 
lies in the future light cone, it follows that both vectors p; and ps» lie in the future 
light cone. We now write 


AE 0 rå Vad 


I, > Å dp, dr, 0 (-pdD 0 (— på) 009 Op) EP ERE 3 


U 
(æ,t) 


; CAR2DD) 
202 
Pi P2 6 (4) X3y 1.7 Be. XTgy 7 Ty FE 
< re g J Y 11/8 ERR å 
<( ) VA 0 9 4 Pi|J, 9 V/ SN 
Pee UENS 
Ai GE pD (A. 21å) 
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If we here introduce "polar coordinates” according to 


Pxo = Tx COSh 0, 


Pr DE SIDE ORE 


we can perform the integrations over the "angles” 6, and obtain 


æ REDER ON 8 
i rs) 0(A ENES 
I, > SAV r, dr; dm” å) 6) mf (- Ze rzer, Ju (-2x E ED 


As the second term in (A. 15) can be obtained from the first one if we replace x,, by 
-T4z7, we have the following result for the integral I in (33): 


Lt Es: 
1=5N FA SSOBSE TEE Fa 0 er m) 4 (20 rer å 
CA224) 
T3y t 505 
NESS N SES WESSS NES ER 
with 
rer) 0(4) 
BEEN EET ea A.24a 
(7973) = dr | Bl. (A. 24a) 
1 1: 
A=—1+4=(17+r12)— — (12 —r2]'. (A. 24b) 
2 16 
Eqs. (A. 24) are identical with Eq. (36a) in the main text. 
Appendix III 
Consider the following square root: 
VED VØEPRNGESE), CAR2SY 


where ft, and t, are two fixed complex numbers, both in the upper half-plane, and 
tis a complex variable. To give this root a well-defined meaning, we introduce two 
cuts in the complex t-plane, one between t; and ft, and the other between —t, and 
—1,, according to Fig. 2. We further require that the square root approaches £? for 
large values of | + |. The imaginary and real parts of this root change their signs on 
the cuts and on the curves given by 


BEER) 0; (A. 26) 


ty 
INS Q=+t0 54 
Qg=-00 NE > Q=—0 
DDIS N SS 
N N 5%, 
N Q OSS 
X LISE 
5 »SS55 
OD t, ”ZAGSK 
Dk SÆLGE + 
00 OSSE SÅN ONNGZÆZEEME 


SD / DS OSS ØST TT 
S tD - 1 Q Fa 
N 


9=-, 


The domain where Re VP(tD>O /s shaded DDS 


The domain where ImVP(t) >O /s shaded me 
Fig. 2. The cuts in the t-plane used to define |/P, (t) in (A. 25). 


where 9 is a real number. When 9 is positive, the imaginary part of VP, (t) changes 
its sign, while the real part changes its sign for negative values of 9. The curves (A. 26) 
pass through the points +t, and +, for eQ = 0. Further, they approach the real and 
imaginary axes asymptotically for large positive values of 9. For large negative values 
of 2, they instead approach the two lines through the origin with the directions 1 + i. 
The general behaviour of these curves is indicated in Fig. 2. In the same figure, we 
have also introduced different shadings for the domains where the real and imaginary 


parts of VP; (t) are positive. 


In the second term in (41) in the main text, VP, (t) is used along a path of inte- 
gration between the points t, and 1, (cf. Figs. 1 and 2) and is defined to have a positive 


imaginary part there. According to Fig. 2, this coincides with the definition of VP, (d) 
given here. 

The first term in (41) contains vær (t)-on a path of integration between ti; 
and t, and with the definition Im VÆR: (DEO ecordimnsetorbissæRe VP, (Gris 
positive on this path, and we conclude 


VP, OD =FiY BID! CAPS 


This is one of the relations used in getting from Eq. (41) to Eq. (42) in the main text. 


INrF9 


t 
N ER 
D= & NG Q $£. 00. /f Q=-—0 
SN > 
N 5 
N ISS 
N KS 
DØ 
% S MUS GS 
7 N SMS OSS 
NNE VS 2 5 
DISSE 4. FrWzao SEGS - 
GG SS DVER RE Mi: EEG 
0 SES SÅN CD i 
SDS NG lg 
D SS N DD 
SSR N 
OSSE N 
SKP N 
0 
ØK N 
g=-07 AN g=-0 
Eg Q=+0 V 


The domain where Re VR) > 0 is shaded PDS 


The domain where ImVB(t) >0 is shaded MAX 
Fig. 3. The cuts in the t-plane used to define VP, (t). 
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Å similar discussion for the polynomial P, (t) = (?? — 13) (? — 12) leads to the cuts 


and curves shown in Fig. 3. From this diagram we find the relation 


V/=P(t)=4i]/P,(t) 


(ATDB8) 


along the path of integration between i, and t,, as well as JÆD: (D—— VP (t) on the 
path between ft; and t,. Eq. (ÅA. 28) is the other relation needed to obtain Eq. (42) 


in the main text. 
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